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ABSTRACT 

We present ALMA and VLA observations of the molecular and ionized gas at 0.1-0. 3" reso- 
lution in the Class protostellar system IRAS 16293-2422. These data clarify the origins of 
the protostellar outflows from the deeply embedded sources in this complex region. Source 
A2 is confirmed to be at the origin of the well known large scale north-east-south- west flow. 
The most recent VLA observations reveal a new ejection from that protostar, demonstrating 
that it drives an episodic jet. The central compact part of the other known large scale flow in 
the system, oriented roughly east- west, is well delineated by the CO(6-5) emission imaged 
with ALMA and is confirmed to be driven from within component A. Finally, a one-sided 
blueshifted bubble-like outflow structure is detected here for the first time from source B to 
the north-west of the system. Its very short dynamical timescale (~ 200 yr), low velocity, and 
moderate collimation support the idea that source B is the youngest object in the system, and 
possibly one of the youngest protostars known. 

Key words: stars: formation — ISM: jets and outflows — stars: individual (IRAS 16293- 
2422) — techniques: interferometric — submillimeter 



1 INTRODUCTION 

Outflows constitute one of the most spectacular, emblematic, and 
relevant early manifestations of the star formation process (e.g. 
Lada 1985; Arce et al. 2007). By removing angular momentum, 
they allow accretion to proceed (Bodenheimer 1995), and through 
their injection of mechanical energy, they can affect the overall evo- 
lution of star-forming regions (Reipurth et al. 1997). Their mor- 
phologies can also provide clues on the properties of the driving 
sources. For instance, the degree of collimation of outflows is be- 
lieved to correlate with the age and mass of the driving protostars 
(Machida et al. 2008; Arce et al. 2007), while multipolar outflows 
are often an indication that the central engine is a multiple system 
(e.g. Carrasco-Gonzalez et al. 2008). 

IRAS 16293-2422 (Ceccarelli et al. 2000; Crimier et al. 2010) 
is located at a distance of 120 pc (Loinard et al. 2008) in the eastern 
part of the Ophiuchus complex. Soon after its discovery (Walker 
et al. 1986), it was found to power two bipolar outflow structures 
(Mizuno et al. 1990). One is oriented in the north-east-south- west 
direction (at a position angle of about 55° and hereafter called the 



NE-SW flow), while the other is oriented nearly east-west (RA. ~ 
110°, hereafter the EW flow). The existence of at least two pro- 
tostellar sources was also established early thanks to high resolu- 
tion radio observations (Wootten 1989, Mundy et al. 1992). These 
objects are known as component A to the south-east, and compo- 
nent B to the north-west, and their projected separation is about 5" 
(Figure 1). While component B remains single even at the highest 
angular resolution available (~ 0.05 arcsec; Rodriguez et al. 2005; 
Chandler et al. 2005), component A breaks up into two centimeter 
sub-components called Al and A2 at resolutions better than about 
0.2" (Figure 1). Because the relative orientation of the A1/A2 pair 
at the time of its discovery was very similar to the direction of the 
NE-SW flow, Al was initially believed to be an ejecta from A2. 
However, analysis of the relative motion of Al and A2 favors a 
scenario where these two sources trace the two stars in a binary sys- 
tem (Loinard 2002; Chandler et al. 2005; Loinard et al. 2007, Pech 
et al. 2010). Sub-arcsecond submillimeter continuum observations 
obtained with the SMA (Submillimeter Array) revealed additional 
structure in component A (Chandler et al. 2005). While the emis- 
sion is dominated by a bright compact source (called Aa) centered 



2 Laurent Loinard et al. 



almost exactly between Al and A2, an additional source (Ab) was 
detected to the north-east of component A (Figure 1). This source 
has not been detected at any other wavelengths so far, and its exact 
nature remains mysterious. 

The NE-SW outflow has long been known to originate in com- 
ponent A (e.g. Hirano et al. 2001; Castets et al. 2001). This associ- 
ation was unambiguously confirmed, and was further constrained, 
when Loinard et al. (2007) witnessed, at radio wavelengths, the 
ejection from A2 of a bipolar system of ejecta which, subsequently, 
moved away from A2 roughly along the direction of the NE-SW 
flow (Pech et al. 2010). The EW flow was recently shown to also 
originate from within component A by Yeh et al. (2008) who pre- 
sented SMA CO(2-l) and CO(3-2) observations at an angular reso- 
lution of a few arcseconds. These SMA data, however, have insuf- 
ficient angular resolution to distinguish between Al and A2 as the 
driving source for the EW flow. Interestingly, there is little evidence 
for outflow activity from component B. J0rgensen et al. (201 1) have 
presented submillimeter (SMA) observations of IRAS 16293-2422 
that resolve the A and B components and do not show strong indi- 
cations for high velocity gas toward B. Yeh et al. (2008) reported 
the detection of a compact blue- shifted CO structure to the south- 
east of source B (their structure b2), but their data was insufficient 
to decide whether or not this was a compact flow from source B. 
Rao et al. (2009) also detected this structure, and argued that it was 
part of an additional outflow driven by source A and oriented at a 
position angle of about 145°. In this Letter, we present combined 
ALMA (Atacama Large Millimeter/submillimeter Array) and VLA 
(Karl G. Jansky Very Large Array) observations of IRAS 16293- 
2422 that help elucidate the origin and nature of the outflows in this 
complex young stellar system. 



2 OBSERVATIONS 
2.1 ALMA data 

IRAS 16293-2422 was observed at A = 0.45 mm with fifteen 12-m 
antennas of ALMA on April 2012, during the science verification 
data program. The 105 independent baselines ranged in projected 
length from 26 to 403 m. The primary beam of ALMA at 0.45 mm 
has a FWHM of 8", so the observations were made in mosaicking 
mode with half-power point spacing between field centers to cover 
both source A and source B. The digital correlator was configured 
in 4 spectral windows of 1875 MHz and 3840 channels each. This 
provides a channel spacing of 0.488 MHz (~ 0.2 km s _1 ) per chan- 
nel, but the spectral resolution is a factor of two lower (0.4 km s _1 ) 
due to online Hanning smoothing. Observations of Juno provided 
the absolute scale for the flux density calibration while observa- 
tions of the quasars J1625-254 and NRAO530 (with flux densi- 
ties of 0.4 and 0.6 Jy, respectively) provide the gain calibration. 
The quasars 3C279 and J 1924-292 were used for the bandpass cal- 
ibration. The data were calibrated, imaged, and analyzed using the 
Common Astronomy Software Applications (CASA) and KARMA 
software (Gooch 1996). The continuum emission as well as many 
spectral lines were detected in these ALMA data (see J0rgensen 
et al. 2012). Here we concentrate on the analysis of the contin- 
uum emission and of the CO(J=6-5; z/=0) at a rest frequency of 
691.47308 GHz. The r.m.s. noise, measured in line-free channels 
on either side of the CO line, is about 50 mJy beam -1 at an angular 
resolution of / . / 32 x O'/IS; -69°. Only a small number of line-free 
channels were averaged to calculate the continuum, so the r.m.s. 
noise of the continuum image is only moderately better (20 mJy 
beam -1 ) than that in the individual spectral channels. 
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Figure 1. Sub-millimeter continuum image of the IRAS 16293-2422 sys- 
tem from ALMA. The contours run from 0.2 to 4 Jy beam -1 by steps 
of 0.2 Jy beam -1 ; the synthesized beam (0 / . / 32 x 0'.'18; -69°) is shown 
at the bottom-right. The noise level is 0.02 mJy beam -1 . The two main 
components (A and B) are labelled, and the direction of the two outflows 
driven from component A are indicated (from Mizuno et al. 1990). The sub- 
millimeter peaks Aa and Ab from Chandler et al. (2005) and the centimeter 
sources Al and A2 are shown. 



2.2 VLA data 

IRAS 16293-2422 was observed with the VLA in its most extended 
(A) configuration on 2011, June 5 and 8 at A = 7 mm, and on 
2011, August 12 at A = 4 cm. In all cases, the absolute flux cal- 
ibrator was 3C286 (J133 1+3030), while the amplitude and phase 
gains were monitored using observations of Jl 625-2527. At 7 mm, 
sixteen contiguous spectral windows, each containing 64 spectral 
channels 2 MHz- wide were recorded simultaneously to produce a 
total bandwidth of 2 GHz (between 40 and 42 GHz). At 4 cm, eight 
contiguous spectral windows, each containing 64 spectral channels 
2 MHz- wide were recorded simultaneously toproduce a total band- 
width of 1 GHz (between 6.8 and 7.8 GHz) The data calibration 
was performed using the CASA software. The data were manu- 
ally inspected and flagged, and calibrated following standard pro- 
cedures. The calibrated visibilities at 4 cm were imaged using a 
nearly uniform weighting scheme (Robust parameter set to -2 in 
CASA) yielding a synthesized beam of 0"44 x 0"16; -11° and a 
noise level of 16 /xJy beam -1 . The visibilities at 7 mm were imaged 
both with natural weighting (Robust = +2) to obtain the highest 
sensitivity (at the detriment of angular resolution), and with inter- 
mediate weighting (Robust = 0) to obtain higher angular resolution 
(at the detriment of sensitivity). The synthesized beam was 0^14 
x 0^12; +69° and the sensitivity was 19 /xJy beam -1 in the former 
case, while they were 0"08 x 0"07; -35° and 31 /iJy beam -1 in 
the latter. While the sensitivity at 4 cm is near the theoretically ex- 
pected value given the integration time, the sensitivity at 7 mm is 
about 50% worse than expected. Additional analysis will be needed 
to understand the origin of this degradation, but we note that in spite 
of this effect, the maps presented here are roughly 5 times deeper 
that the best previous 7 mm images (Rodriguez et al. 2005). 



1 A second spectral band, corresponding to the frequency range from 4.2 
to 5.2 GHz was observed simultaneously with the 4 cm data presented here, 
but will not be discussed in this Letter. 
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Figure 2. VLA continuum images of IRAS 16293-2422. (a) 7 mm image 
reconstructed with natural weighting to optimize sensitivity. The synthe- 
sized beam (0'.'14 x 0'/12; +69°) is shown at the top-right of the panel, and 
the sensitivity is 19 /Jy beam -1 . The contours are at 70, 140, 300, 600, 900, 
etc. /Jy beam -1 . For reference, the position of the sub-millimeter peaks Aa 
and Ab are shown as magenta crosses. Note the faint extension to the NE co- 
incident with Ab and the extension of the sub-millimeter ALMA image, and 
the faint emission associated with the ejecta A2a and A2/3. (b) 7 mm image 
reconstructed with intermediate weighting. The synthesized beam (O'.'OS x 
/ . / 07; -35°) is shown at the top-right of the panel, and the sensitivity is 31 
/Jy beam -1 . The contours are at 150, 300, 600, 900, etc. /Jy beam -1 . Note 
the new ejecta A2S. (c) 4 cm image where all the sources and ejecta are 
seen. The synthesized beam (0"44 x r !l6; -11°) is shown at the top-right 
of the panel, and the sensitivity is 16 /Jy beam -1 . The first contour and the 
contour spacing are at /Jy beam -1 . The double arrow shows the direction 
(at P. A. 70°) along which the ejecta A2cn, A2/3, and A2S were launched. 



3 RESULTS AND DISCUSSION 

3.1 Continuum emission and the NE-SW flow 

The sub-millimeter continuum emission (A = 0.45 mm) mea- 
sured by ALMA reveals the usual double source structure of 
IRAS 16293-2422 (Figure [TJ. Component B is resolved but com- 
pact, with a peak brightness temperature of about 160 K. The total 
flux measured here with ALMA (~ 12.2 Jy) is consistent with the 
trend observed at lower frequencies that the spectral index of source 
B is between 2 and 2.5 from centimeter to submillimeter wave- 
lengths (Chandler et al. 2005). Component A is elongated along the 
NE-SW direction, with its peak at the position of source Aa, and the 
faint extension to the NE corresponding to the position of source 
Ab reported by Chandler et al. (2005). Interestingly, the centimeter 
sub-components Al and A2 are not identifiable in this image, al- 



though the angular resolution would be sufficient to separate them. 
The peak brightness temperature of the sub-millimeter emission in 
component A is about 30 K. The overall properties of the ALMA 
image of component A suggest that the sub-millimeter continuum 
is dominated by the circumbinary envelope rather than by the in- 
dividual disks around the protostars. This is not particularly sur- 
prising since the individual circumstellar disks are expected to be 
heavily tidally truncated in this very compact binary system (pro- 
jected separation of order 40 AU). In contrast, the continuum at 7 
mm is dominated by free-free radiation from the bases of the ther- 
mal jets and enables us to pin-point the embedded sources Al and 
A2 (Figure [2}. This demonstrates the complementarity of ALMA 
and the VLA to study young stellar sources. 

When it is reconstructed to optimize sensitivity, the 7 mm 
emission from source A also reveals a faint elongation toward 
the NE similar to that seen in the sub-millimeter map, and ex- 
tended emission surrounding the sources Al and A2. Like the sub- 
millimeter continuum, this extended emission likely traces the cir- 
cumbinary envelope. The emission associated with the extension of 
source A to the north-east (i.e. with source Ab from Chandler et al. 
2005) contributes roughly 1.4 mJy at 41 GHz, 0.5 Jy at 305 GHz 
(Chandler et al. 2005) and 3 Jy at 696 GHz (all values with relative 
uncertainties of order 50%). This yields a spectral index of 2.7 ± 
0.5, consistent with that of the more compact emission from source 
A (Chandler et al. 2005). We argue that this elongation traces faint 
extended emission from the circumbinary envelope (possibly re- 
lated to dust accumulated in that direction by the NE-SW flow), 
rather than a new protostellar source in the system. 

Both the 7 mm and the 4 cm continuum maps (Figure [2]) re- 
veal a new compact source to the NE of source A2. This is strik- 
ingly similar to the situation encountered in the 2006 image of 
IRAS 16293-2422 at A = 1.3 cm (Loinard et al. 2007) when a bipo- 
lar ejection from A2 was seen for the first time. The bipolar ejection 
that occurred in 2006 produced the sources A2a and A2(3 that are 
still clearly seen in the 4 cm image (Figure [2]:) and, faintly, in the 7 
mm map (Figured). They are moving away from A2 at projected 
velocities of several tens of km s _1 (Pech et al. 2010). We interpret 
the new source to the NE of A2 (which we shall call A25) as a new 
ejection from A2. In the high resolution 7 mm image, there is a faint 
extension to the SW of A2 which might correspond to the bipolar 
counterpart of A25. We note that in the first image following the 
ejection of the A2a/A2/3 pair (Loinard et al. 2007), the SW ejecta 
(Ala) was not immediately detached from A2 either. 

IRAS 16293-2422 has been regularly monitored at centime- 
ter wavelengths since the late 1980s (Chandler et al. 2005) and 
the 2006 ejection was the first to be detected. The new ejection 
reported here suggests that source A2 might be entering a period 
of enhanced outflow activity. Regular observations in the coming 
years and decades will be required to confirm this possibility, and 
will enable a detailed study of the kinematics and behavior of the 
ejecta. Interestingly, A2a, A2/3, A2S, and A2 itself are well aligned 
(Figured) but the position angle (~ 70°) of the line joining them is 
somewhat larger than the position angle (~ 55°) of the large-scale 
outflow known to originate from A2 (e.g. Figure 5 in Mizuno et al. 
1990). This suggests that the jet driven by source A2 is precessing 
as a consequence of the A1-A2 binarity. We note, finally, that there 
is no detectable CO(6-5) emission associated with the NE-SW flow 
in the ALMA observations (Figure©. This is in agreement with the 
results of Yeh et al. (2008) who also failed to detect such emission 
in their SMA CO observations at arcsecond resolution. 
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Figure 3. Integrated intensity of the weighted velocity (moment 1) color 
map of the CO(6-5) emission from source A overlaid in contours with the 
0.45 mm continuum emission (black thick line) and the integrated intensity 
line emission (moment 0) of the CO(6-5) (grey thin line). The black con- 
tours are the same as in Figure [T] (0.2 to 1.2 Jy beam -1 in steps of 0.2 Jy 
beam -1 ). The grey contours are from 10% to 90% with steps of 10% of the 
peak of the CO line emission; the peak is 3.2 x 10 4 Jy Beam -1 km s -1 . The 
color-scale bar on the right indicate the LSR velocities in km s -1 . The syn- 
thesized beam of the continuum image is shown in the bottom left corner of 
the image. 



3.2 The EW flow 

Strong CO(6-5) is detected in the ALMA data around the position 
of source A. In agreement with the results of Yeh et al. (2008), it is 
largely concentrated to an EW structure with blueshifted emission 
to the east and redshifted emission to the west; This emission traces 
the central part of the well-known EW (PA ~ 110°) flow (Figure[3]). 
In spite of the high angular resolution of the ALMA observations, 
the origin of the flow stall cannot be traced back to a specific source 
within component A. 



3.3 A compact outflow from source B 

As reported by Yeh et al. (2008) and confirmed by Rao et al. (2009), 
there is very strong blue-shifted CO emission to the south of source 
B (the location of component b2 in Yeh et al. 2008). Figure HI dis- 
plays the first moment CO(6-5) map of this component (in colors) 
overlaid with the 0.45 mm continuum image (in contours). It is 
clear that the CO emission defines a 3 /r -long bubble-like structure 
oriented along the south-east-north- west direction (at RA. ~ 130°) 
with source B at its north-west apex. This structure points from 
source B toward source A, but there appears to be no material con- 
nection (i.e. no bridge of emission) between the two. In particular, 
its south-east apex (corresponding to the point nearest to A) is lo- 
cated about 2" to the north-west of source A. We conclude that this 
structure is associated with source B, and unrelated to source A. 
We should point out, here, that there is a significant amount of ex- 
tended CO emission around the systemic velocity of IRAS 16293- 
2422 (vi sr ~ 4 km s _1 ) in the ALMA observations. This emission is 
poorly recovered by the interferometer, and its structure cannot be 
assessed; we filtered it out by ignoring the velocity channels around 
4 km s" 1 , but note that our fluxes (as given, for instance, in Figure 
|4]) are clearly underestimated since they do not include the ambient 
gas. Interestingly, Rao et al. (2009) recently reported on the possi- 
ble detection of a new outflow in IRAS 16293-2422 seen in SiO(8- 
7) along a position angle very similar to that of the CO structure 
seen here. The structure seen in SiO, however, is most prominent 
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Figure 4. Integrated intensity of the weighted velocity (moment 1) color 
map of the CO(6-5) emission from source B overlaid in contours with the 
0.45 mm continuum emission (black thick line) and the velocity scale of 
the CO(6-5) (grey thin line). The black contours are the same as in Fig- 
ure [T] (starting at and in step of 0.2 Jy beam -1 ). The color-scale bar on the 
right indicate the LSR velocities in km s _1 . The synthesized beam of the 
continuum image is shown in the upper left corner of the image. 



south of source A, and entirely at positive velocities. While there 
is SiO(8-7) emission in the general direction of the CO structure 
reported here, it is at +4 km s _1 rather than at -4 km s _1 . Higher 
resolution SiO observations will clearly be needed to establish the 
relation (if there is any) between the SiO structure reported by Rao 
et al. (2009) and the CO structure reported here. 

To our knowledge, ours is the first direct indication for the ex- 
istence of an outflow driven by source B. Interestingly, there is no 
strong redshifted counterpart to the north-west. While this might re- 
flect an intrinsic asymmetry of the flow, it could also at least partly 
reflect the structure of the circumstellar environment. As discussed 
by Pineda et al. (2012) and Zapata et al. (in prep.), at sub-millimeter 
wavelengths, both the CO lines and the continuum emission toward 
source B are optically thick, so the red-shifted emission might be 
at least partly hidden from view. Alternatively, IRAS 16293-2422 
might be located close to the back side of the cloud where it is lo- 
cated. Besides explaining the lack of a red- shifted lobe, this would 
be consistent with the very high optical depth toward source B. 

Since the systemic velocity of source B is about 3 km s _1 (J0r- 
gensen et al. 2011) while the most negative velocities reached by 
the outflow are ~ -8 km s" 1 , the (radial) expansion velocity of the 
bubble relative to source B is about 10 km s" 1 . At the distance of 
Ophiuchus, the angular extent of the bubble (3") corresponds to 5.4 
x 10 15 cm (360 AU), and the dynamical age of the flow is of or- 
der 200 years. While this calculation might only provide a crude 
estimate of the true age of the structure, similar calculations for 
the other outflows in the system (Mizuno et al. 1990) yield much 
larger dynamical ages, of order 10 4 years. This is consistent with 
the commonly held view that source B is the youngest object in 
IRAS 16293-2422 (e.g. Chandler et al. 2005). It could be argued 
that the true angular extent of the outflow is significantly larger than 
reported here, because the CO(6-5) line might only pick the most 
highly excited gas in the vicinity of source B. This is unlikely be- 
cause lower J CO observations (e.g. Yeh et al. 2008) did not reveal 
any conspicuous emission on larger scales along the direction of 
this structure. 

The 7 mm continuum emission around source B has usually 
been interpreted as a nearly face-on disk (Rodriguez et al. 2005). 
Recent ALMA spectroscopic observations further revealed unam- 
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biguous evidence for infall and rotation (Pineda et al. 2012, Zapata 
et al. in prep.) toward that source, and the rotation is also consistent 
with a nearly face-on orientation. In this configuration, the flow 
driven by source B ought to be oriented nearly along the line of 
sight, and the expansion velocity would be mostly radial. This sug- 
gests that the true expansion velocity of the flow is not significantly 
larger than the observed 10 km s _1 . 

The overall morphology of the outflow driven by source B is 
qualitatively different from that of typical flows powered by low- 
mass young stars. In particular, there is no evidence for a collimated 
jet-like structure along the symmetry axis. Instead, the highest (i.e. 
most negative) velocities are associated with the inner layer of the 
structure, while the outer layers are at less negative velocities (typi- 
cally at about -2 km s _1 ; see Figure |4}. This behavior suggests that 
source B drives a moderately collimated wind, that impinges into 
the stationary circumstellar medium, in a situation resembling that 
described theoretically by Shu et al. (1991). 

In summary, the flow driven by source B has a number of pe- 
culiar properties. It is fairly slow (~ 10 km s _1 ), outburst-like, and 
only moderately collimated. Those are characteristics theoretically 
expected from flow driven by adiabatic (first) cores (Machida et al. 
2008, Price et al. 2012; see also Chen et al. 2010, Enoch et al. 2010, 
Dunham et al. 2011, Pineda et al. 2011 for observational aspects). 
However, the evidence for strong accretion reported by Pineda et 
al. (2012) and Zapata et al. (in prep.) and the fact that component 
B is such a bright submillimeter continuum source would be more 
consistent with the idea that it has already entered the protostel- 
lar phase. In any case, the properties of the outflow driven by this 
object do suggest that it might be one of the youngest known pro- 
tostars. 



4 CONCLUSIONS AND PERSPECTIVES 

In this Letter, we used combined ALMA and VLA observations to 
examine the outflow activity in the very young (Class 0) protostel- 
lar system IRAS 16293-2422. The well-known NE-SW and E-W 
large-scale molecular flows are confirmed to be driven from within 
component A to the south-east of the system. The NE-SW flow 
can unambiguously be associated with source A2 which might be 
entering a phase of enhanced activity. The ALMA CO(6-5) obser- 
vations have enabled us to discover a new outflow toward source B 
with peculiar properties: it is highly asymmetric, resembles an out- 
burst, is fairly slow (10 km s _1 ), and lacks a jet-like feature along 
its symmetry axis. In addition, its dynamical age is only about 200 
years. We argue that these properties result from the extreme youth 
of source B, which might be one of the youngest known proto- 
stars. Given the short dynamical age of the structure, significant 
morphological changes are expected to take place over the next 
few decades. Observing this time evolution would help constrain 
the formation and early evolution of molecular outflows in gen- 
eral, and would have a fundamental impact on our understanding 
of these objects. 
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